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ABSTRACT 

We study high redshift structure formation and reionization in a ACDM universe under the assumption 
that the spectral power index of primordial density fluctuations is a function of length scale. We adopt 
a particular formulation of the "running" spectral index (RSI) model as suggested by the combined 
analysis of the recent WMAP data and two other large-scale structure observations. We carry out high 
resolution cosmological simulations and use them to study the formation of primordial gas clouds where 
the first stars are likely to form. 

While early structure forms hierarchically in the RSI model, quite similarly to the standard power-law 
ACDM model, the reduced power on small scales causes a considerable delay in the formation epoch 
of low mass (~ lO 6 Af ) "mini-halos" compared to the ACDM model. The abundance of primordial 
star-forming gas clouds in such halos also differs by more than an order of magnitude at z > 15 between 
the two models. The extremely small number of gas clouds in the RSI model indicates that reionization 
is initiated later than z < 15, generally resulting in a smaller total Thomson optical depth than in the 
ACDM model. 

By carrying out radiative transfer calculations, we also study reionization by stellar populations formed 
in galaxies. We show that, in order to reionize the Universe by z ~ 7, the escape fraction of ultraviolet 
photons from galaxies in the RSI model must be as high as 0.6 throughout the redshift range 5 < z < 18 
for a stellar population similar to that of the local Universe. Even with a top-heavy intial mass function 
representing an early population of massive stars and/or an extraordinarily high photon emission rate 
from galaxies, the total optical depth can only be as large as r e ~ 0.1 for reasonable models of early 
star-formation. The RSI model is thus in conflict with the large Thomson optical depth inferred by the 
WMAP satellite. 

Subject headings: cosmology: theory - early Universe - reionization of the Universe 



1. INTRODUCTION 

The origin of matter density fluctuations in the Universe 
is one of the fundamental problems in cosmology. The 
so-called standard theory of structure formation posits 
that the present-day clumpy appearance of the Universe 
developed through gravitational amplification of an ini- 
tially very smooth, but perturbed matter distribution. Re- 
cent observations provide a consistent picture of the large- 
scale matter distribution in the early Universe; fluctua- 
tions arise from adiabatic perturbations whose statistics 
are described by a Gaussian field, as predicted by popular 
inflationary theories. The inflationary models also pre- 
dict a scale-invariant power spectrum of primordial den- 
sity fluctuations. In fact, the generic inflationary predic- 
tion is nearly but not precisely scale-invariant density fluc- 
tuations; i.e. the fluctuation power spectrum scales as 
P(k) (x k n with n « 1 (see Lyth & Riotto 1999 for a 
recent review). The tilt, or the exact value for n, has at- 
tracted considerable attention in theories of the physics of 
the early Universe (Covi & Lyth 1999), in astronomical 
observations (Croft et al. 2002; Hannestad et al. 2002; 
Seljak, McDonald & Makarov 2003) and in the context of 



structure formation (Zentner & Bullock 2002). 

The spectral index of the primordial power, n, may even 
be a function of length scale. This "running" of the spec- 
tral index modulates the relative amplitudes of density 
fluctuations on large to small scales. Interestingly, the 
combined analysis of the first-year Wilkinson Microwave 
Anisotropy Probe (WMAP) data, the 2dF galaxy red- 
shift survey, and Lyman-a forest observations favors a 
cosmological model with a "running" of the primordial 
power spectrum (Spergel et al. 2003; Peiris et al. 2003). 
Although the WMAP data alone gives a best-fit power 
law ACDM model with the spectral index n ~ 1 and 
the fluctuation amplitude erg = 0.9, the Lyman-a forests 
observations (Croft et al. 1998, 2002; McDonald et al. 
2000; Gncdin & Hamilton 2002) consistently favor a low- 
amplitude model on scales 1 < k [Mpc -1 ] < 10, which can 
be reconciled if the primordial power spectrum has a mild 
tilt (n < 1) or a negative running (dn/dhik < 0). For 
a given normalization on a large scale, models with nega- 
tive running predict progressively reduced linear power on 
small scales (Kosowsky & Turner 1995; Hannestad et al. 
2002). 

Suppressing linear power on small, galactic to sub- 
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galactic scales may offer a resolution to the often-claimed 
problems of Cold Dark Matter (CDM) models (e.g. Moore 
et al. 1999). It could alleviate apparent discrepancies 
between the observed structure of dark halos and predic- 
tions from CDM models (Zentner & Bullock 2003). On 
the other hand, too much reduction of small-scale power 
may result in a conflict with other observations. Dalai & 
Kochanek (2002) and Chiba (2002) argue that the abun- 
dance of dark matter substructure in several galactic lens 
systems inferred from observed flux anomalies is consis- 
tent with predictions from the standard ACDM model 
(see however Bullock & Zentner 2002 for an attempt to 
constrain some variant models). In models with reduced 
small scale power, nonlinear objects with characteristic 
mass 10 6 — 10 9 M Q form late, and hence reionization due 
to photons from stellar sources formed in low mass sys- 
tems is expected to occur late (Somcrville, Bullock & Livio 
2003). To the contrary, the measurement of TE polariza- 
tion by WMAP suggests a large Thomson optical depth 
r e = 0.17 ± 0.04, implying that reionization took place as 
early as z ro i on ~ 17 (Kogut et al. 2003). 

Somerville et al. (2003) argue that models like the RSI 
one we consider here are inconsistent with the inferred 
early reionization epoch. Yoshida et al. (2003b) studied 
early structure formation in a warm dark matter model 
and showed that suppressing the linear power even on very 
small scales (k > 100 Mpc -1 ) can considerably affect the 
formation of high redshift star- forming gas clouds, making 
early reionization by massive, metal-free stars unlikely to 
occur early on. Chiu, Fan & Ostrikcr (2003) argue that a 
large optical depth r e > 0.1 may require exotic radiation 
sources other than a normal stellar population. Avclino 
& Liddlc (2003) also conclude that running of the primor- 
dial power systematically delays the reionization epoch. 
Overall, these theoretical studies suggest that, if the large 
Thomson optical depth observed by the WMAP satellite 
is confirmed, the standard power-law (n « 1) ACDM 
model is favored over other models with reduced small 
scale power. Thus, the WMAP results appear to suggest 
a somewhat contradictory cosmological model. Since the 
major reason to favor power-law ACDM in this context 
is the inferred early epoch of reionization, it is clearly im- 
portant to study and compare quantitatively details of the 
reionization process in the two cosmologies, the RSI model 
and a power law ACDM universe. 

In the present paper, we study early structure for- 
mation in RSI and ACDM universes using cosmologi- 
cal Af-body/hydrodynamic simulations. By combining 
these dynamical models with radiative transfer calcula- 
tions, we compute the epoch of reionization by various 
stellar sources in the two cosmologies and compare the 
Thomson optical depth to that measured by the WMAP 
satellite. 

In particular, at high redshifts, z £ 15, we consider the 
implications of an early population of massive, metal-free 
stars (the "first stars"), which have been hypothesized to 
form in "mini-halos" (~ 10 6 M Q ). We refer to the sites 
where these stars can form as "primordial gas clouds;" i.e. 
dense regions of cold gas which accumulate through the 
action of cooling by molecular hydrogen. In order to iden- 
tify these primordial gas clouds, our simulations follow the 
non-equilibrium chemical reactions of a chemically pristine 
gas and molecular hydrogen cooling. In an earlier study 



(Yoshida et al. 2003a), we found that non-equilibrium ef- 
fects associated with dynamical heating play a critical role 
in regulating the supply of molecular gas, contrary to the 
assumptions commonly made in semi-analytic studies of 
this process. Hence, it is necessary to follow the dynamics 
of structure formation explicitly to reliably estimate the 
abundance of early stars. It is also important to specify 
the source location accurately when carrying out radia- 
tive transfer caculations. Unlike semi-analytic methods, 
direct hydrodynamic simulations enable us to locate plau- 
sible star-forming sites robustly in a cosmological volume. 
As we discuss later in the present paper, the distribution 
of the sources strongly affects the overall topology of the 
ionized regions. 

We show that the reduced small scale power in the 
RSI model causes a considerable delay in the formation 
epoch of low mass halos and the primordial gas clouds 
within them compared to a ACDM counterpart. In view 
of this result, it appears unlikely that massive stars in 
mini-halos can contribute significantly to reionization in 
cosmologies with a running spectral index, at least for 
those with parameters similar to what is inferred from the 
joint WMAPext + 2dFGRS + Lyman-a forest analysis 
(Spergel et al. 2003). To study the impact of ordinary 
( "Population II" ) stars formed in galaxies and protogalax- 
ies, we carry out ray tracing calculations of reionization by 
stellar sources. For an initial mass function (IMF) similar 
to that in local galaxies, reionization can be completed at 
sufficiently early times to account for the Lyman-a optical 
depths measured in quasars at z > 6, but in the case of the 
RSI model only if a large fraction of ionizing photons are 
able to escape from galaxies. Furthermore, because reion- 
ization occurs late in this context, the Thomson optical 
depth is too small to be compatible with that suggested 
by the WMAP measurements. 

Finally, we consider the possibility that the production 
rate of ionizing photons could be boosted in galaxies and 
protogalaxics with a top-heavy IMF. In this case, Thomson 
optical depths approaching that favored by the WMAP 
analysis are possible, but only if the parameters describ- 
ing these stars are pushed to extreme values. 

The paper is organized as follows. We first describe the 
cosmological simulations we use to study the abundance 
of high-redshift star-forming systems in §2. We present 
the simulation results in §3. In §4, we discuss the possible 
effect of the finite size of our simulation volume. In §5, 
the abundance of primordial gas clouds and implications 
are described. There, we also give a brief description of 
the basic relation between reionization epoch and the to- 
tal Thomson optical depth. In §6 we show the results of 
radiative transfer calculations in which we consider ordi- 
nary stellar populations in galaxies and protogalaxies as 
radiation sources. A summary and discussion are given in 
§7- 
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Fig. 1. — The linear power spectra for the power law ACDM 
model (solid line) and for the RSI model (dashed line). In the bot- 
tom panel we plot the variance A 2 (A;) linearly extrapolated to the 
present epoch for the two models. Note the difference in the plotted 
range in k between the top and the bottom panels. The arrow in 
the bottom panel indicates the length scales our simulations probe. 



2. THE iV-BODY /HYDRO DYNAMIC SIMULATIONS 

We use the parallel Tree-PM solver GADGET2 com- 
bined with smoothed particle hydrodynamics (SPH) in our 
simulations. Our treatment of SPH employs the conserva- 
tive entropy formulation of Springel & Hcrnquist (2002), 
which offers several distinct advantages over previous ver- 
sions of SPH. In particular, because the energy equation 
is written with the entropy as the independent thermo- 
dynamic variable, as opposed to the thermal energy, the 
'p&V term is not evaluated explicitly, reducing noise from 
smoothed estimates of e.g. the density. By including terms 
involving derivatives of the density with respect to the par- 
ticle smoothing lengths, this approach explicitly conserves 
entropy (in regions without shocks), even when smooth- 
ing lengths evolve adaptively, avoiding the problems noted 
by e.g. Hcrnquist (1993). Furthermore, this formulation 



moderates the overcooling problem present in earlier for- 
mulations of SPH (Yoshida et al. 2002; see also Pearce et 
al. 1999, Croft et al. 2001). 

The simulations follow the non-equilibrium reactions of 
nine chemical species (e~, H, H + , He, He + , He ++ , H2, 
H^, H~) using the reaction coefficients compiled by Abel 
et al. (1997). We use the cooling rate of Galli & Palla 
(1998) for molecular hydrogen cooling. Further simula- 
tion details are found in Yoshida et al. (2003a, b). The 
simulations employ 2 x 324 3 particles in a comoving vol- 
ume of 1 Mpc on a side. We work with A-dominated cos- 
mologies with matter density f2 m = 0.3, cosmological con- 
stant f^A = 0.7 and the Hubble constant at the present 
time h — 0.7 in units of 100km s _1 Mpc~ 1 . We set the 
baryon density to f2b = 0.04 and the normalization pa- 
rameter to as = 0.9. The transfer functions are computed 
by the Boltzmann code of Sugiyama (1995) for the adopted 
cosmology and initial conditions are realized as described 
in Yoshida, Sugiyama & Hcrnquist (2003). The primor- 
dial power spectrum for the ACDM case is assumed to 
be P(k) oc k, whereas that of the RSI model is given by 
P(k) oc k n » with 



ns(fco) + ^dlnT ln V/ i(1 



(1) 



Specifically, we choose /c =0.05 Mpc 1 , n s (fc )=0.93, and 
dn s /d\nk=-Qm, as indicated by the WMAPext + 2dF- 
GRS + Lyman-a analysis (Spergel et al. 2003). Note that 
we adopt the same normalization, erg = 0.9, for both the 
power law ACDM model and the Running Spectral Index 
model. Hereafter, we refer to the conventional power law 
ACDM model simply as the "ACDM model," and to the 
other as the "RSI model" . 

Figure ^ shows the linear power spectra P(k) for the 
two models. In the bottom panel, we plot the variance 
per unit logarithmic interval in k, A 2 (fc) = fc 3 P(fc)/(27r 2 ), 
to better describe the difference between the two models 
in the range 0.1 < fcfMpc- 1 ] < 1000. Since the input 
power spectra are normalized at a length scale larger than 
the simulation boxsize, the amplitudes of the density fluc- 
tuations realized in the RSI simulation are smaller than 
in the ACDM simulation on all length scales. Hence, the 
density fluctuations on large scales (k ~ 10) in the RSI 
model should evolve similarly to those in ACDM but with 
a slight delay in time. We note that, while the boxsize of 
our simulations is sufficiently large to ensure that the fun- 
damental fluctuation mode in the simulations remains in 
the linear regime at the early epochs we simulate, it may 
not sample a fair volume of the Universe. We will discuss 
finite boxsize effects in §4. 

Figure [21 shows the projected gas distribution at z = 22 
(top panels) and at z = 17 (bottom panels) for the two 
models. We choose these output redshifts such that the 
density fluctuation on the largest length scale in the RSI 
model at z = 17 resembles that of the ACDM model at 
z = 22. As expected, the RSI model appears as a delayed 
ACDM model on these scales. Also, rich small-scale struc- 
ture develops early on only in the ACDM model, as can 
be inferred from the linear power spectrum (Figure 
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Fig. 2. — The projected gas density distribution at 2 = 22 (top 
panels) and z = 17 (bottom panels) for the ACDM (left) and RSI 
(right) models. 




3. HALO ABUNDANCE 

Since the typical masses of the halos that host early 
primordial gas clouds are 5 x 10 5 - 10 7 M© (Tegmark et 
al 1997; Machacek et al. 2001; Yoshida et al. 2003a), 
the abundance of such small mass halos is an important 
quantity to compare. In order to quantify the difference 
between the two models, we measure the mass function 
of dark matter halos. Whereas locating dark halos in 
pure iV-body simulations is relatively simple using con- 
ventional techniques such as the friends-of-friends (FOF) 
algorithm, there are a few subtleties in identifying and as- 
signing properties to halos in A^-body/hydrodynamic sim- 
ulations. Since our main objective is to compare the halo 
abundance between the two models, we avoid ambiguities 
by adopting a simple procedure, as follows. 

We locate halos by running a FOF groupfinder with link- 
ing parameter b = 0.2 on the dark matter particles in our 
simulations. We then assign the total mass to each group 
as the number of linked member particles multiplied by the 
scaled particle mass m par t = md m + m gas , where nidm is 
the dark matter particle mass and m gas is the gas particle 
mass. We compute the Press-Schechter (PS) mass function 
to which we input the linear power spectrum for the dark 
matter component. (Note that, for our simulations, sepa- 
rate transfer functions are used for the dark matter and gas 
components in setting-up the initial conditions [for details, 
see Yoshida, Sugiyama & Hcrnquist 2003].) We compare 
the mass function of the halos identified in this manner at 
z = 17 and z = 22 in Figure For both models, the 
measured mass function is overall well fitted by the PS 
mass function (the solid line for the ACDM model and the 
dashed line for the RSI model). Again, the RSI model at 
z = 17 resembles the ACDM model at z = 22 with the 
abundance of low-mass halos being appreciably small. 

*Only in Figure 131 and Figure ITT1 we use units that include the Hubble 
mass and (h~ 1 Mpc) — 3 for abundance. 



Fig. 3. — The mass function for the ACDM model (left) and for 
the RSI model (right) at 2 = 17 and 2 = 22. The solid line is the 
Press-Schechter mass function for the ACDM model, and the dashed 
line is that for the RSI model. 

The halo abundance in the RSI model appears to devi- 
ate slightly from the PS mass function. Although a simi- 
lar feature is also seen in the ACDM model, the deviation 
from the PS mass function is more clearly noticeable for 
the RSI model. We discuss this issue together with finite 
boxsize effects in the next section. 



4. MASS VARIANCE 

It is well known that structure in CDM models grows 
hierarchically; i.e., smaller mass objects form earlier and 
merge to form bigger objects. While this can be easily 
inferred from the variance A 2 (fc) for the ACDM model, 
which has larger amplitudes on smaller length scales (see 
Figure nj, the situation is less clear for the RSI model. 
As Figure shows, the variance in the RSI model ap- 
pears approximately constant over a range of scales 10 < 
fcpVTpc -1 ] < 1000, indicating that nonlinear objects with 
widely different masses could form nearly at the same 
epoch. It is also worth mentioning that it is not clear 
whether the Press-Schechter theory can be used for such 
models. For the comparison shown in Figure |3| we cal- 
culated the PS mass function by formally inputting the 
linear power spectrum for the RSI model. Although the 
agreement shown in Figure is reasonably good, it is not 
trivially expected. The flatness of the variance of the RSI 
model deserves particular attention also for another rea- 
son; the hierarchical nature of the formation of dark mat- 
ter halos has a strong influence on the formation of primor- 
dial gas clouds (Yoshida et al. 2003a). Hence it is impor- 
tant to verify whether or not small mass (~ 10 6 M Q ) non- 
linear objects in the RSI model grow hierarchically early 
on. For this purpose, it is more appropriate to work with 
a properly defined mass variance rather than the power 
spectrum. 

constant h, to be consistent with those commonly used; h~ 1 Mq for 
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Fig. 4. — The mass variance computed from the power spec- 
tra linearly extrapolated to the present epoch. The thick solid line 
is for the ACDM model, and the thick dashed line is for the RSI 
model. The thin lines show the effective mass variance as defined 
in equation for each model. 

We compute the mass variance as 

a 2 (M) = J^[ P(k)W 2 (kR)k 2 dk, (2) 
2n z J 

where the top-hat window function is given by W{x) = 
3(sin(a;)/a; 3 — cos(x)/x 2 ). Figure 0] shows er(M) for the 
two models. In the plotted range, the mass variance in the 
RSI model monotonically increases toward smaller mass 
scales, indicating that less massive objects form at ear- 
lier epochs; i.e., structure formation is still expected to 
be "bottom-up" . Note that the RSI model adopted in 
the present paper predicts an effective slope of the power 
spectrum smaller than -3 for k ^> 200. On very small 
length scales (k > 1000), or equivalently on low mass 
scales [M < 1000Af Q ), structure formation should then 
proceed differently from the "bottom-up" picture. Owing 
to the flattening of the variance, the characteristic mass 
scale of the first non-linear objects may not be well de- 
fined. Halos with masses 10 — 1000M© should collapse 
first and larger objects will be formed through mergers of 
these small halos. Simulations with a substantially higher 
mass resolution are needed to reveal how halos with such 
very low masses form, although it will not be important 
for the formation of halos with masses ~ 10 6 M Q that are 
crucial for primordial gas cloud formation. 

Let us now examine the consequences of the finite box 
size of our simulations. We quantify the effect by comput- 
ing an effective variance 

a 2 s (M) = ^J k P{k)W 2 {kR)k 2 dk, (3) 

where we set the lower boundary k m i n to exclude the con- 
tribution from large-scale density fluctuations which are 
missing in our simulations. We simply set fc nml = 27r/Lbox 
with Lbox being the simulation box sidelength. The thin 
lines in Figure 0] are the effective mass variance computed 
in this manner. For both the models, a c g(M) is 
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Fig. 5. — The number of halos with mass greater than 7 X 10 5 Mq 
per cubic mega-parsec (comoving) volume. The solid line is the an- 
alytic estimate obtained from the Press-Schechter mass function for 
the ACDM model, and the dashed line is for the RSI model. 

slightly smaller than a(M) on the relevant mass scales 
(~ 10 5 — 10 6 Af Q ), and thus a systematic delay in the char- 
acteristic formation epoch (or a slightly lower abundance 
at a given epoch) is expected. This may account for the 
small discrepancy in the mass function plotted in Figure 
|3J Since we use linearly extrapolated power spectra to 
compute a(M) and <r e ff(M), the precise effect of a finite 
boxsize remains still uncertain. The apparent difference 
shown in Figure 0] should be interpreted as qualitatively 
representing the finite boxsize effect. Note that the RSI 
model is more affected by the finite boxsize because the 
relative fractional power on large scales is bigger than in 
the ACDM model (see Figure^. Systematic studies using 
larger box simulations are clearly needed to address this 
issue further. 

By constructing a halo merger history using 30 simu- 
lation outputs dumped frequently between 17 < z < 25, 
we also verified that almost all the halos that have masses 
greater than 7 x 10 5 M Q at z = 17 grew monotonically 
in mass by repeated mergers and continuous accretion of 
smaller mass elements, quite similarly to what is found 
in the ACDM model at somewhat higher redshifts. As 
long as we are concerned with the abundance of low-mass 
halos, the results of our small box simulations appear to 
be in reasonable agreement with the Press-Schechter mass 
function. In Figure [S] we show the number of halos with 
masses greater than M crrt = 7 x 10° Mq. We chose the 
threshold mass M cr i t based on the results of Yoshida et 
al. (2003a) for star-forming systems. Although the in- 
complete sampling of the halo mass function owing to the 
finite boxsize is appreciable at high redshift, z > 27 for 
the ACDM model and z > 20 for the RSI model, overall 
the halo abundance is in agreement with the PS predic- 
tion. Thus, we conclude that the statistics we discuss in 
the present paper, such as the abundance of primordial gas 
clouds, are not significantly affected by the finite simula- 
tion boxsize. The apparent discrepancy in Figure [S] may 
be partly due to some inaccuracies of the analytic mass 
function itself (see Jenkins et al. 2001; Reed et al. 2003), 
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although the simulation results of Jang-Condell & Hern- 
quist (2001) indicate that the PS mass function is a good 
approximation on these mass scales. An important point 
to note is that both the analytic estimate and the sim- 
ulation results suggest that the number of of halos with 
mass ~ 10 6 M©, which are crucial for the primordial gas 
cloud formation, is an order of magnitude smaller in the 
RSI model. 

Finally, we mention that a promising way of compen- 
sating for the missing large-scale power would be to em- 
ploy the mode-adding procedure developed by Tormen & 
Bertschinger (1996) and by Cole (1997). The method can 
effectively extend the dynamic range of cosmological simu- 
lations, and thus will be useful for studies of early structure 
formation in a cosmological context. 

5. CONTRIBUTION OF THE "FIRST STARS" TO 
REIONIZATION 

Numerical simulations (e.g. Abel, Bryan & Norman 
2002; Bromm, Coppi & Larson 2002) suggest that the first 
stars in the Universe were unusually massive and formed 
in "mini-halos" at redshifts z > 20, as a consequence of 
cooling by molecular hydrogen. It is customary to re- 
fer to early, massive stars as Population III. In what fol- 
lows, we will distinguish between massive stars which may 
have formed in mini-halos by molecular cooling (the "first 
stars") from those which could have originated in larger 
halos where gas can cool by atomic processes. While the 
nature of the "first stars" has been studied in some detail, 
the stellar populations expected in slightly more massive 
halos (the "first galaxies") is unclear. With modest un- 
certainty, we can quantify the importance of the former 
to reionization, but numerous assumptions are required to 
account for "non-standard" IMFs in early protogalaxies, 
as we discuss later in the paper. 

5.1. Primordial gas cloud formation 

Owing to the delayed formation of "mini-halos" in the 
RSI model, the formation of primordial gas clouds is also 
expected to occur late. It is in these regions that sig- 
nificant quantities of molecular gas can accumulate, lead- 
ing to the formation of massive, metal-free stars of the 
type modeled by e.g. Abel, Bryan & Norman (2002) and 
Bromm, Coppi & Larson (2002). We emphasize that non- 
equilibrium processes not typically accounted for in semi- 
analytic estimates, like the dynamical heating effect iden- 
tified by Yoshida et al. (2003a), significantly influence the 
abundance of these star-forming regions. 

In FigureElwe plot the number of gas clouds identified in 
the two simulations. We define groups of cold (T < 500K), 
dense (jih > 500cm~ 3 ) gas particles as "gas clouds". We 
locate the gas clouds by running a FOF groupfinder to 
the simulation gas particles with a small linking parame- 
ter b = 0.05. We then discard from the groups gas par- 
ticles that do not satisfy the above conditions. Figure 
clearly shows that the total number of gas clouds in the 
simulated volume differs by more than an order of mag- 
nitude between the two models in the redshift range plot- 
ted. At z — 17, we identified 66 gas clouds in the CDM 
model, whereas there is only one gas cloud found in the 
RSI model. We continued the RSI run down to z = 15 
to see if there is a rapid increase in the number of gas 
clouds. Although the increase in the number of gas clouds 



at z ~ 14 — 17 appears large, the total number is still much 
smaller than that of the ACDM model. As we argued in 
section 4, some halos grow very rapidly in the RSI model, 
which causes considerable perturbation to the formation 
and the growth of the primordial gas clouds within them. 
This may partly explain the extremely small number of 
gas clouds at z ~ 15 in the RSI model. It is interesting 
that the number of the gas clouds in the RSI model is 
even smaller than, or at most comparable to, the warm 
dark matter case we studied in Yoshida et al. (2003b). 
For the specific model 
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Fig. 6. — The number of primordial "gas clouds" where massive, 
metal free stars can form by molecular hydrogen cooling. 

employed in Yoshida et al. (2003b) that assumes a single 
massive (3OOM0), metal-free star is formed per gas cloud, 
reionization by the "first stars" is not completed at z > 15 
in the RSI model. In fact, with the very small number 
of sources, the ionized volume fraction rises only up to a 
few percent by z — 15. This is clearly in disagreement 
with the WMAP result of z IC i on ^17. Since the estimated 
uncertainty in the WMAP result of the reionization epoch 
is rather large (Kogut et al. 2003), it still remains to 
be seen whether the particular RSI model we consider is 
compatible with the WMAP data. Nevertheless, it is clear 
that the running of the power spectrum is not favored in 
terms of early reionization. 

5.2. Reionization epoch and the Thomson optical depth 

The differential optical depth to Thomson scattering in 
a small redshift interval dz is 

dr e = cr T n e (z) c*^- dz, (4) 

where ctt = 6.65 x 10~ 25 cm 2 is the Thomson scattering 
cross-section, n e (z) is the mean electron number density 
at z, and c is the speed of light. The mean electron number 
density is given by 

n e (z) = n e (0)(l + z) 3 x e (z), (5) 

where x e is the ionization fraction. For a flat ACDM 
model, H(z) = H a [n A + ft m (l + zf] 1 ' 2 . 
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Fig. 7. — The optical depth to Thomson scattering to z r cion when 
abrupt and complete reionization at z r cion is assumed. 

Putting x e = 1 for an interval between z\ and z%, we 
obtain the contribution as 



2 <tt n e (0)c 

3(0^)1/2 



+ + 



f2/ 



1/2 



(6) 



Consider a simple case where abrupt and complete 
reionization occurs at z. Using the WMAP results for the 
cosmological parameters, f2b^ 2 =0.0224 and 17 m /i 2 =0.135, 
and assuming a plasma of primordial abundance with, for 
example, fully ionized hydrogen and singly ionized helium, 
we obtain the total optical depth to redshift z: 



dr P w 0.0023 



([1 + z] 3 + 2.7) 



1/2 



1.93 



(7) 



Figure d shows the total optical depth against reioniza- 
tion epoch at z = -z r cion- The WMAP result is indicated 
at r e = 0.17 (solid line) with a 1-cr range 0.17 ±0.04 (dot- 
ted line). The actual reionization history could be more 
complicated, but Figure provides the useful insight that 
a large Thomson optical depth, say r e > 0.13, can be 
achieved only if reionization takes place at an early epoch, 



> 14. 



Our simulation of early structure formation in the RSI 
model showed that there are only 5 gas clouds at z = 15 
in the simulated volume. As Yoshida et al. (2003b) 
demonstrate, approximately 100 very massive stars must 
be turned on within ~1 recombination time in a (1 Mpc) 3 
volume to cause complete reionization. Owing to the 
strong radiative feedback from the first stars, the num- 
ber of stars formed in the primordial gas clouds is likely to 
be limited to one (Omukai & Nishi 1999; Oh et al. 2001). 
Hence the number of "first stars" in the RSI model is far 
too small at z > 15 to reionize the Universe. The epoch of 
complete reionization is certainly much later than z ~ 15, 
implying a small total optical depth r e <C 0.13. 

In summary, primordial gas cloud formation which relies 
on molecular hydrogen cooling in low-mass, "mini-halos" 
appears to be extremely inefficient at z > 15 in the RSI 
model and thus other sources, such as early protogalaxies, 



must contribute to reionization for this model to be com- 
patible with the WMAP results. Furthermore, it is impor- 
tant to establish whether or not ordinary stellar sources 
in galaxies in the RSI model can account for reionization 
early enough to even be consistent with measurements of 
the Lyman-a optical depth in quasars at z > 6. 

6. REIONIZATION BY STARS IN GALAXIES 

In the previous section, we have shown that low-mass 
mini-halos and primordial gas clouds which can harbor 
massive, metal free stars (the "first stars") do not form in 
sufficient numbers at high redshift in the RSI model to ac- 
count for reionization at z ^ 15. As we emphasize in §5.2 
and Figure 7, reionization must occur early if the WMAP 
measurement of the Thomson optical depth is accurate. 

The number of larger mass (~ 10 7 — 10 9 M©) halos is 
also predicted to be small for the RSI model. Thus, reion- 
ization by stars in systems that are the likely progenitors 
of galaxies will be less efficient than for a ACDM model, 
for similar assumptions about the stellar populations. In 
this section, we study reionization by stars in galaxies and 
proto-galaxies with total masses larger than ~ 4 x 10 7 M Q , 
where the gas can cool rapidly by atomic line transitions. 
The nature of the stars forming early in these environ- 
ments is uncertain. It is possible, for example, that the 
IMF will be top-heavy, at least initially, perhaps increas- 
ing the production rate of ionizing photons per unit stellar 
mass compared to a conventional IMF inferred for stellar 
populations in the local Universe. 

In what follows, we begin by considering star formation 
in halos in the RSI model by assuming that the stellar 
population is similar to that locally and is unevolving. 
Namely, we adopt a Salpeter (1955) IMF and normalize 
the star formation rate (SFR) so that it reproduces the 
Kennicutt (1998) law. For simplicity, we refer to this de- 
scription as a "Population II" SFR. In §6.4 we examine the 
overall effect of relaxing these assumptions. However, it is 
already clear that the RSI model will have difficulty ac- 
counting for the WMAP estimate of the Thomson for any 
reasonable assumptions about star formation simply be- 
cause structure forms too late in this scenario to produce 
a large r e , as suggested by Figure 7. 

To study how reionization proceeds in the RSI model, we 
run multi-source radiative transfer simulations using the 
technique of Sokasian et al. (2001, 2003a, b). In order to 
locate radiation sources and to compute the IGM density 
field, we utilize existing outputs of a large volume ACDM 
simulation from Springel & Hernquist (2003b), normaliz- 
ing the star formation rate expected for the RSI model 
using the analytic description of Hernquist & Springel 
(2003), which can be applied to any cosmology. 

6.1. The star-formation rate in the RSI model 

First, we compute the "generic" star-formation rate for 
the RSI model using the analytic model of Hernquist & 
Springel (2003). Their approach describes the formation 
of dark matter halos using an analytic halo mass function 
(Sheth & Tormen 1999; 2002), and employs simplified pre- 
scriptions for the cooling of gas and star-formation which 
are commonly used in semi-analytic models of galaxy for- 
mation (Kauffmann et al. 1999; Springel et al. 2001; 
Yoshida et al. 2002). This analytic formalism has been 
tested in detail by comparing its predictions to those from 
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Fig. 8. — The star-formation rate density in galaxies ("Popula- 
tion II") as a function of redshift for the ACDM (solid line) and RSI 
(dashed line) models. 

cosmological simulations that include hydrodynamics and 
star formation (Springel & Hernquist 2003b). A thorough 
description of the model is given by Hernquist & Springel 
(2003) and references therein. 

In particular, Hernquist & Springel (2003) show that 
under rather general circumstances the star-formation rate 
density is well-approximated by the simple form 



where 



p»(z) = p*(0) 



xO) 



X" 



1 + a( X - 1)« exp(/3 x r ) ' 



H(z) 
H 



2/3 



(8) 



(9) 



In equation JHJl , the factors (3 and r are fixed by the under- 
lying cosmological power spectrum, the choices for p and 
q reflect the interplay between radiative cooling in halos 
and the expansion of the Universe, and the normalizing 
constants p*(0) and a depend on the rate at which cold, 
dense gas is converted into stars. 

The form of equation (jHJ) is motivated by two compet- 
ing effects which regulate star formation. At sufficiently 
high redshifts, collapsed gas is dense and cooling times are 
short so that star formation is only limited by the grav- 
itational growth of massive halos. At low redshifts, the 
cooling time in halos increases owing to the decline in the 
mean density of the Universe by cosmic expansion. In this 
regime, the supply of star-forming gas is throttled by the 
expansion of the Universe. 

Based on the multiphase description of star-forming gas 
described by Springel & Hernquist (2003a), Hernquist & 
Springel (2003) show that the appropriate values for the 
parameters in equation JSJ) for the ACDM model employed 
here are p = 2, q = 3, r = 7/4, a = 0.012, (3 = 0.041, and 
/5,(0) = O.O13Af /yr/Mpc 3 . 

For the RSI cosmology, we use the RSI linear power 
spectrum as input and follow the procedure described in 
Hernquist & Springel (2003). We find that the resulting 
SFR density is again well-fitted by equation © for the 



choices p = 1.5, q = 3.05, r = 1.85, a = 0.01, (3 = 0.057, 
and p*(0) = O.O176Af /yr/Mpc 3 . 

In Figure El we compare the SFR density estimated in 
this manner for the two cosmologies. The RSI model pre- 
dicts a significantly smaller SFR than the ACDM model, 
in qualitative agreement with the semi-analytic estimates 
of Somerville et al. (2003). We use the forms for the SFR 
density in our ray-tracing calculatiosn only within the re- 
stricted redshift range 5 < z < 20 to compute radiation 
source luminosities. 

6.2. Source luminosities 

The second step is to prescribe the luminosity of the 
sources ("galaxies") located in the simulations. Specif- 
ically, we use the outputs from Run Q5 of Springel & 
Hernquist (2003b), which simulates a cube of 14.3 Mpc on 
a side. We first calibrate the SFRs of individual sources 
identified in the simulation, and then convert these into 
luminosities. As in Sokasian et al. (2003a), we use the de- 
rived SFR (equation |HJ as a reference to match the total 
SFR of the sources. We carry out the calibration proce- 
dure after selecting sources in the simulation outputs in 
the following two ways. 

6.2.1. Dim source model 

In this approach, we identify all the halos found in each 
output as "galaxies" and assign luminosities to them based 
on the actual star-formation rates which are computed in 
a dynamically consistent way in the cosmological simula- 
tion. The minimum mass of the sources is set by the reso- 
lution of the simulation, which is M m i n = 4.3 x 10 7 /i -1 Mq 
in our case. The SFRs of the sources with masses close 
to the resolution limit are corrected in the same manner 
as in Sokasian et al. (2003a, §4.2.1). We then scale the 
SFRs of all the galaxies by multiplying by a constant fac- 
tor Cl(z) such that the total SFR within the simulated 
volume matches that for the RSI model obtained from 
equation (JSJ. Note that the scaling factor Cl(z) is kept 
constant for all the galaxies within one simulation output, 
but it is time-dependent. 

6.2.2. Massive halo model 

Since the halo abundance itself is lower in the RSI model 
than in the ACDM cosmology, it may be more appropri- 
ate to populate a smaller number of objects as radiation 
sources in the RSI model. Motivated by this, we first com- 
pute the number density of halos with virial temperature 
greater than 10 4 K in which the gas can cool by atomic 
line transitions. The threshold mass is given by 



M min = 10 s 



1 



11 



-3/2 



Mr. 



(10) 



Then we use the halo mass function to compute the ex- 
pected abundance of halos with mass greater than M m j n . 
Let us denote this number by N somcc (z) . In the simulation 
output produced at z, we select N souxce (z) halos in order 
of their mass. Luminosities are assigned to the selected 
sources in proportion to their mass, L oc Afhaio such that 
the total SFR summed over the selected sources matches 
the global SFR given by equation © , with the parameter 
values appropriate for the RSI model. 
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6.2.3. Escape fraction 

Finally, the ionizing flux from the sources is modulated 
by the "escape fraction," / csc , which describes the proba- 
bility that individual photons can escape from the galaxy 
in which they were produced into the surrounding inter- 
galactic medium (IGM). Sokasian et al. (2003a) studied 
reionization by stellar sources in the ACDM universe by 
adopting several values for / esc . Based on a detailed com- 
parison with the observational results of the z — 6.28 
quasar (SDSS 1030.10+0524) in Becker et al. (2001), 
they conclude that / csc = 0.2 is most consistent with 
the observed Lyman-a flux transmittance. We take this 
/csc = 0.2 model in Sokasian et al. (2003a) as our fiducial 
case for the ACDM model. For the RSI model, our zeroth 
order requirement is that reionization should be completed 
by z = 7. After several experiments, we found that set- 
ting / csc > 0.6 is necessary for the RSI model to satisfy this 
condition. We thus set / esc = 0.6 as our basic choice for 
the RSI model and compare the results with the fiducial 
ACDM model. 

6.3. Numerical results for "Population II" star formation 

We compute the total volume filling factor of the ionized 
medium Qhii from the results of our ray-tracing calcula- 
tions. Figure shows the resulting ionized volume frac- 
tion, and Figure ITU1 shows the Thomson optical depth r e 
as a function of redshift. In the figures, we show only the 
results for our "massive-halo model" of source selection 
for the RSI model. Comparing the results for the dim- 
source model and the massive-halo model, we found that 
there is no substantial difference between the two in ei- 
ther T e (z) or Qhii(z). We note, however, that the overall 
morphology of the ionized regions looks somewhat differ- 
ent between the two cases, because there are more sources 
in the dim-source model that reside in less dense regions. 
Nevertheless, for global quantities such as the total optical 
depth, the two models are found to be quite similar. 

On the other hand, the results for the ACDM and the 
RSI models differ substantially. As Figure shows, the 
ionized volume fraction in the RSI model remains small 
until z rsj 9, and then rapidly increases close to unity at 
z ~ 7. At later epochs, z <J 9, the abundance of high 
mass (~ 10 9 M Q ) halos rapidly increases, and the total 
star-formation (and hence the photon production) is dom- 
inated by such large halos. This can also be understood in 
terms of the halo mass function at z = 7 shown in Figure 
1111 By this time, the mass functions for the two models are 
already quite similar, with the difference being only a con- 
stant factor of roughly two (compare with the much larger 
difference at z = 17 at smaller mass scales shown in Fig- 
ure^. Namely, the RSI model quickly catches up with the 
ACDM model in the abundance of halos around z ~ 7—11 
(see also the increasing SFR in Figure|SJ), which causes the 
rapid increase in the ionized volume fraction as seen in Fig- 
ure [3] Reionization is nearly completed (Qhii > 0.95) at 
z = 8 in the ACDM model, whereas it occurs at z ~ 7.3 
in the RSI model. However, as shown in Figure the 
resulting total optical depth r e = 0.055 for the RSI model 
is still small, owing to the small Qhii at earlier epochs. 
We emphasize that a large Thomson optical depth can be 
achieved only if reionization occurs early (§6). In other 
words, after reionization is completed, the additional 
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Fig. 9. — The ionized volume fraction for the ACDM model with 
/esc = 0.2 (solid line) and for the RSI model with / eS c = 0.6 (dashed 
line), for "Population II" stars. 
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of our simulations for the ACDM model with f CBC = 0.2 (solid line) 
and for the RSI model with / eS c = 0.6 (dashed line) for "Population 
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photon supply does not matter in the total optical depth 
as long as the IGM is kept almost fully ionized thereafter, 
as implied by Figure 7. (Note that here we discuss only 
hydrogen reionization.) 

Our fiducial ACDM case with / csc = 0.2 predicts a larger 
total optical depth than the RSI model with / csc = 0.6 
by Ar e ~ 0.01. The difference mostly comes from the 
fact that Qhii is still quite small at z > 9 in the RSI 
model. Interestingly, none of the "conventional" scenarios 
we have considered for "Population II" star formation re- 
produces the claimed high optical depth from WMAP of 
0.13 <r e < 0.21 (1-ct; Kogut et al. 2003). It is clear that 
these models are inconsistent with the WMAP measure- 
ment and thus may need, for example, additional radiation 
sources to ionize the IGM early on. 

In related work, we have shown that an early generation 
of massive, metal-free stars in mini-halos supplemented by 
"ordinary" star formation at lower redshifts to account for 
the spectra of high redshift quasars can yield a Thomson 
optical depth approaching that of the WMAP measure- 
ment in a ACDM cosmology (Yoshida et al. 2003b, 



10 



CD 
O 



10" 



10 £ 



10' 



10 1 - 



10" 



10" 



10" 



10" 



I I ' 


i i 

LCDM _ 




RSI 


z=7 





10' 



10 c 



10" 



10 10 10 11 10 12 



Msun/h 



7. 



Fig. 11. — The Press-Schechter mass function at z - 

Sokasian et al. 2003b), along the lines suggested by 
Haiman & Holder (2003) and Cen (2003). However, it 
does not appear that a similar conclusion will obtain for 
the RSI model because of the low abundance of mini-halos 
at high redshifts, as we found in the case of a warm dark 
matter model with particle mass 10 keV (Yoshida et al. 
2003b). 

6.4. Star formation in galaxies with boosted photon 
emission 

The results described above indicate that the RSI model 
is incompatible with a high Thomson optical depth like 
that suggested by the WMAP measurements. A remain- 
ing loophole in this argument is that the stellar population 
in galaxies at early times may have differed from that in 
our underlying model so that a larger number of ionizing 
photons were produced from galaxies at z > 10 than we as- 
sumed above (e.g. Ciardi, Ferrara & White 2003; Wyithe 
& Loeb 2003; Sokasian et al. 2003a). We first discuss 
such possibilities quantitatively, and then show numerical 
results for a model that enhances the photon production 
rate in galaxies. 

For a Salpeter-type power-law IMF, the hydrogen ioniz- 
ing photon emission rate of zero-metallicity stars is larger 
than that of a Population II (i.e. having metals) coun- 
terpart only by about 50% (Thumlinson & Shull 2000; see 
Leitherer et al. 1999 for low-metallicity cases.). The stellar 
IMF in the first galaxies may have had even an extreme 
top-heavy shape, with massive stars exclusively forming 
out of a primeval gas (Schwarzschild & Spitzer 1953; Mat- 
suda, Sato & Takeda 1969; Yoshii & Saio 1986; Larson 
1998; Abel, Bryan & Norman 2002; Bromm, Coppi & Lar- 
son 2002). Very massive Population III stars are efficient 
emitters of UV photons, with the hydrogen ionizing pho- 
ton emission rate per stellar mass being up to 10-20 times 
higher than for ordinary stellar populations (Bromm, Ku- 
dritzki & Loeb 2001; Schaerer 2002). It is thus conceivable 
that, in early galaxies, the photon emission rate per stellar 
mass was ~ 10 times higher than that of the "Population 
II" case we considered above. 

More specifically, the predicted star-formation rate for 
the RSI model shown in Figure [SJ is 
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Fig. 12. — The Thomson optical depth for the RSI model in which 
the photon emission rate is "boosted" according to equation 1111 . 
The filled square is the WMAP result with its l-cr error bar. 

~ 5 x lO~ 4 M0yr -1 Mpc -3 at z — 15. Assuming a con- 
stant conversion factor of 10 53 s _1 ionizing photons per 
star- formation rate of 1 MQyr^ 1 for an ordinary stel- 
lar population with solar metallicity (Madau, Haardt & 
Rees 1999), we obtain a total rate of ~ 5 x 10 49 ionizing 
photons per second per cubic comoving 1 Mpc volume 
at z = 15 for our fiducial RSI model. If all the formed 
stars have a mass of 300M Q , the star-formation rate of 
~ 5 x 10 _4 M©yr^ 1 Mpc -3 corresponds to the formation of 
5 such stars every three million years within the volume. 
The calculation of Bromm et al. (2001) suggests that 
a total emission rate of 2.5 x 10 51 ionizing photons per 
second, a 50 times higher photon emission rate, can be 
maintained in this case, assuming such very massive stars 
are continuously formed. (Note that these massive stars 
are typically short-lived, with a main-sequence lifetime of 
about two million years.) We thus conclude that a 100 
times higher photon production rate than assumed in our 
original RSI model serves as a conservative upper limit. 

We examine the overall effect of such an additional pho- 
ton supply by "boosting" the photon emission rate of the 
sources in the RSI massive-halo model. Sokasian et al. 
(2003a) found that boosting the photon emission rate at 
high redshifts progressively shifts the reionization epoch 
(see their Figure 13). Motivated by their result, we em- 
ploy a simple prescription which increases the emission 
rate of individual sources according to 



N, 



ph.eff. 



exp(0.53 (z - 6)) x N ] 



ph j 



(11) 



where iV p h and iV p h, e ff. & re the original and the 
boosted photon emission rates, respectively. With this 
parametrization, the ionizing photon emission rate of the 
galaxies steeply increases toward higher redshift, modeling 
either a transition of the stellar IMF or a contribution from 
massive Population III stars or both. This rather extreme 
model requires effectively more than a 10 times higher pho- 
ton emission rate at z > 10, and more than 100 times 
higher at z > 15. We carry out radiative transfer simula- 
tions for the RSI model with this luminosity "boost" . Fig- 
ure El shows the resulting r e (z). The reionization epoch 
is found to be shifted to earlier times, z ~ 10, than in 
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the case considered in §6.3, and consequently r e is greatly 
increased compared with the fiducial RSI model. How- 
ever, the total optical depth for this boost model reaches 
T e ~ 0.12 at z = 0, still lower than the 1-a range of the 
WMAP result. 

7. SUMMARY 

Our numerical simulations of early structure formation 
show that primordial gas cloud formation in low-mass 
"mini-halos" is very inefficient at z > 15 in the RSI model, 
making it unlikely that the "first stars" contribute signif- 
icantly to reionization in this scenario. Using radiative 
transfer calculations, we also show that reionization by 
ordinary ("Population II") stellar sources in galaxies is 
completed late in the RSI model if we employ conven- 
tional models of star-formation. In order to be compatible 
with the epoch of the end stage of reionization inferred 
from high redshift quasar observations, the photon escape 
fraction from galaxies must be large, / csc ^0.6. We found 
that the resulting total Thomson optical depth in this case 
is r e « 0.055, in apparent conflict with the recent WMAP 
measurements of CMB polarization. 

In order for the total optical depth to be as large as 
T e > 0.1, reionization must be completed before z ~ 10 
(Figure©. Ciardi et al. (2003) and Wyithe & Loeb (2003) 
suggest that employing an effectively top-heavy IMF may 
be necessary to cause early reionization. Cen (2003) argues 
that a large Thomson optical depth could be obtained if 
the "Population III era" is prolonged so that massive stars 
continuously form until quite low redshift. If a large vol- 
ume of the IGM could be kept chemically pristine until low 
redshift (z ~ 10) due to the very low global star-formation 
rate in the RSI model, the stellar IMF in galaxies may be 
top-heavy (Schneider et al. 2002; Mackey, Bromm & Hern- 
quist 2003). Unfortunately, little is presently known about 
star-formation in a primordial gas in high mass halos -the 
"first galaxies"- and clearly needs to be addressed using 
numerical simulations as well as theoretical modeling (Lar- 
son 1998; Omukai 2001; Oh & Haiman 2002). We have 
explored the possibility of enhancing the photon emission 
rate by employing a simple prescription. The model pa- 
rameters were chosen such that the photon emission rate 
is close to a plausible maximum bound. Our analysis has 
shown that, whereas the extreme boost model does pre- 
dict an earlier reionization epoch z rc i on ~ 10 and a larger 
total optical depth r e ~ 0.12, the result is still marginally 
inconsistent with the WMAP data. 

Recent numerical simulations by Bromm, Yoshida & 
Hernquist (2003) show that early Population III super- 
novae quickly pollute the surrounding IGM with metals. 
If many massive population III stars are formed in "mini- 
halos" and explode as hyper- energetic supernovae, the 
global metal-enrichment of the IGM is quickly achieved 
to such a degree that a transition of the stellar IMF is 
caused at a very early epoch (Yoshida, Bromm & Hern- 
quist 2003). Then, the large photon production rate at 



z > 6 we hypothesized as coming from metal-free stars 
cannot be maintained for a sufficiently long time, and the 
reionization history in the RSI model will be just as shown 
in Figure El with a small total optical depth (r e < 0.06) 
as shown in Figure ITU1 

In the end, we are left with a contradiction: the run- 
ning of the primordial power spectrum favored by the 
WMAPext + 2dF + Lyman-a analysis and the high 
Thomson optical depth measured by the WMAP satel- 
lite appear to be inconsistent with one other. We have 
discussed various possibilities to resolve the conflict, and 
conclude that a rather radical solution or a fine-tuned com- 
bination, if any is possible, of the proposed resolutions are 
necessary. Models with a weaker running of the primor- 
dial power spectrum will alleviate the discrepancy slightly, 
but the overall results would not be changed. If the high 
Thomson optical depth is confirmed and the running of 
the primordial power spectrum turns out to be real, it may 
be necessary to invoke exotic radiation sources other than 
stellar populations or some ionization mechanisms. Such 
scenarios include decaying particles (Sciama 1982; Dodel- 
son & Jubas 1992; Hansen & Haiman 2003) or formation 
of early mini-quasars (Eisenstein & Loeb 1995; Sasaki & 
Umemura 1996). 

The prospects for observationally resolving this dis- 
crepancy in the near future appear bright. Data from 
planned CMB polarization experiments, by the continued 
operation of WMAP, and post-WMAP observatories such 
as Planck will pin down a precise value for r e . Lyman- 
a forest observations exploiting a large sample of SDSS 
quasars will place a tighter constraint on the matter power 
spectrum on large scales, as will the SDSS galaxy redshift 
survey. In the longer term, it may be possible to map 
out the evolution of reionization from redshifted 21 cm 
emission using instruments such as the Square Kilometer 
Arrayt or the Low Frequency Array"'" (e.g. Madau et al. 
1997; Tozzi et al. 2000; Iliev et al. 2002; Ciardi k, Madau 
2003). For example, Furlanetto et al. (2003) show that 
frequency fluctuations can be used to distinguish between 
models with a different number of reionization epochs. 
Analyses of these various high-precision data promise to 
provide a more complete picture of the matter density 
distribution in the early Universe over a wide range of 
scales and its relationship to the formation of stars and 
galaxies. 
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